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We investigate the crossover behavior between special and ordinary surface transitions in three- 
dimensional semi-infinite Ising-like systems with random quenched bulk disorder. We calculate the 
surface crossover critical exponent the critical exponents of the layer, ai, and local specific heats, 
an, by applying the field theoretic approach directly in three spatial dimensions (d = 3) up to the 
two-loop approximation. The numerical estimates of the resulting two-loop series expansions for the 
surface critical exponents are computed by means of Pade and Pade-Borel resummation techniques. 
We find that $, ai, an obtained in the present paper are different from their counterparts of pure 
Ising systems. The obtained results confirm that in a system with random quenched bulk disorder 
the plane boundary is characterized by a new set of critical exponents. 

T3 

C ■ I- INTRODUCTION 

o 

In recent decades the remarkable progress in understanding of the critical behavior of real physical systems was 
achieved from application of the powerful field theoretical methods and renormalization group (RG) approach to 
the analysis of these systems. It was allowed to perform with higher accuracy the numerical analysis of critical 
exponents and universal amplitude combinations for bulk phase transitions [1]. Moreover, these methods have given 
the possibility to investigate the preasymptotic behavior of infinite systems. A series of field-theoretical methods 
developed and tested in the studies of bulk phase transitions, have been extended to study the critical behavior of 
. systems with boundaries. General reviews on surface critical phenomena are given in Refs. [3-5]. 

The presence of surfaces, which are inevitable in real systems, leads to additional complications. A typical model to 
study the critical phenomena in real physical systems restricted by a single planar surface is the semi-infinite model 
[3] . As it is known from field-theoretic analyses of continuum </> 4 model [4] , the influence of the surface is possible to 
take into account by a quadratic surface term with coeficient Co, which describes the enhancement of the interactions 
at the surface, and additional surface fields hi. There are different surface universality classes, defining the critical 
behavior in the vicinity of boundaries, at temperatures close to the bulk critical point (r = (T — T c )/T c — > 0). Each 
bulk universality class, in general, divides into several distinct surface universality classes. Three surface universality 
classes, called ordinary (cq — > oo), special (co = c* p ) and extraordinary (erj — > -co), are relevant for our case [4,5,7]. 
They meet at a multicritical point (to 2 ,, cq) = (toq c , c* p ), which corresponds to the special transition and is called the 
special point [6]. 

Most theoretical studies usually concentrate their attention on the investigation of critical behaviour distinctly at the 
fixed points cq — ±oo and cq — c* p , respectively. At the present time is very well developed theory of critical behavior 
of individual surface universality classes for pure isotropic systems [7-10,5,11] and systems with quenched surface- 
enhancement disorder [12-14]. General irrelevance-relevance criteria of the Harris type for the systems with quenched 
short-range correlated surface-bond disorder were predicted in [12] and confirmed by Monte-Carlo calculations [13,15]. 
Besides, the investigation of the critical behavior of the semi-infinite systems with random quenched bulk disorder 
at the ordinary and the special transitions have been studied by us [16,17]. The obtained results [16,17] have shown 
that such systems are characterized by the new set of surface critical exponents in comparison with the case of pure 
systems. 

As it is known, experimental systems are typically characterized by the parameters different from the fixed point 
values. But, the understanding of the situation in the crossover regions between different transitions is less complete. 
Investigations of the crossover behavior between different surface universality classes for pure isotropic systems have 
been published in a series of papers [18,10,19,20,11,21]. However, at the present time is open the question about the 
picture in the crossover regions between the different transitions for the semi-infinite systems with random quenched 
bulk disorder. In the present paper we restrict our attention to the simplest case with hi = and investigate 
crossover behavior between special and ordinary transitions for semi-infinite Ising-like systems with random quenched 
bulk disorder. Here it should be mentioned that from the whole class of O(N) symmetric TV-vector models in d 
dimensions only Ising model is the one of primary interest, because it satisfies Harris criterion for the specific heat 
exponent a(d) > [22]. 
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The proposed calculations are very important because they allow to understand the phenomenon of adsorption of 
fluid mixtures in contact with a wall, as well as the critical behavior of so-called dilute magnets with the surfaces, 
which can be prepared by mixing an (anti)-ferromagnetic material with nonmagnetic one. 

The calculations are performed by applying field theoretic approach directly in d = 3 dimcntions up to the two-loop 
order approximation. The numerical estimates of the resulting two-loop series expansions for the surface crossover 
exponent $ from the special to the ordinary transition and surface critical exponents of the layer, ai, and local specific 
heats, an, are computed by means of the Pade [23] and Pade-Borcl [24] rcsummation techniques. We find that $, 
ai, an obtained in the present paper are different from their counterparts of pure Ising systems. 

II. MODEL 

The Hamiltonian of the semi-infinite model under consideration with random quenched bulk disorder is given by 

H = ~\ J ijPiPj s i s 3- J 'ij S i S h 

<i,j>£bulk <i.j>£surface 

where Si and Sj are classical m-component spins located at the lattice sites i and j; a nerest-neighbor bond < i,j > 
is said to belong to the surface region if both i £ surface and j <G surface, in other cases they belong to the bulk 
region. The bulk interaction potential Jij have the parallel to the plane translational invariance in the underlying 
lattice. The surface interaction potential will never be invariant with respect to lattice translations parallel to the 
plane or perpendicular to it. The random site variable pi and pj have the probability distribution 

P{pi) =pS(p t - l)+p'5(pi), 

where p = 1 — p is the concentration of nonmagnetic impurities. As it is known, there are two possible ways to 
analyze the above random model. The first way is connected with direct averaging over random disorder using the 
method introduced by Lubensky [25]. The second possibility is to perform the replica trick n — > 0, as it was first done 
in the renormalization group (RG) calculations by Grinstein and Luther [26] . Performing calculation in the spirit of 
the method introduced by Grinstein and Luther it is possible to show that random model (2.1) is thermodynamically 
equvalcnt to the n - vector cubic anisotropic model with effective Hamiltonian of the Landau-Ginzburg- Wilson (LGW) 
type in semi-infinite space at the replica limit n — > 

H{$) =l°°dzj d d -\[\ | V$ | 2 +\ml | $ | 2 

+ ^oE^ + ^o(|^| 2 ) 2 )]+ / d d -\\cj\ (2.2) 

i— 1 

where <p{x) is an n- vector field with the components (fii(x), i — l,...,n. Here m 2 , is the "bare mass" representing 
linear measure of the temperature difference from the critical point value. The values u and v are the usual "bare" 
coupling constants uq < and vq > 0. The constant Co is relates to the surface enhancement, which measures the 
enhancement of the interactions at the surface. It should be mentioned that the d-dimensional spatial integration 
is extended over a half-space Mi_ = {x=(r, z) e M d | r e _/R d_1 ,z > 0} bounded by a plane free surface at z = 0. 
The fields <fri(r,z) satisfy the Dirichlet boundary condition in the case of ordinary transition: (pi(r,z) — at z = 
and the Neumann boundary condition in the case of special transition: d z <pi(r,z) = at z — [7,9]. The model 
defined in (2.2) is restricted to translations parallel to the boundaring surface, z = 0. Thus, only parallel Fourier 
transformations in d — 1 dimensions take place. It should be mentioned that LGW model works good for sufficiently 
low spin dilution 1 — p as long as system is not too close to the percolation limit. 

In order to investigate the critical behavior in the crossover region and to calculate the crossover exponent $ we 
should consider correlation functions with insertions of the surface operator </> 2 

AT M Li 

G^ M ^)({xj,{r,},{R,}) = (n^n^rU^))' (2 - 3) 

i=l j=l 1=1 

which involve N fields </>(xj) at distinct points Xj(l < i < N) off the surface, M fields (f)(rj,z = 0) = ^> s ( r j) at 
distinct surface points with parallel coordinates rj(l < j < M), and Li insertions of the surface operator i(/> 2 (R;) 
; I - / • /.; :. 



The corresponding parallel Fourier transform of the full free propagator takes form 

(2.4) 



G(P,^') = ^ 



e - K0 \z-z'\ _ CQ - KQ_ e - Ko ( z+z >) 



Co + K 



with the standard notation k = y 'p 2 + m 2 . Here, p is the value of parallel momentum associated with d — 1 
translationally invariant directions in the system. 

III. RENORMALIZATION 

The formulation of the renormalization process for the random systems introduced by Grinstcin and Luther [26] 
is essentially the same as in the 'pure' case [4,11]. From other side, as it is known from the theory of semi-infinite 
systems [4,7,8,11], the bulk field </>(x) and the surface field 4> s (r) should be reparameterized by different uv-finitc 

1/2 I/2I/2 

renormalization factors [4,11] Z,p(u,v) and Zi(u,v). Thus we have <fi — <j>R an d 4>s = Z^ (j) s ,R- Besides, 
introducing the additional surface operator insertions i0 2 (R;) requires additional specific renormalization factor Z$i 

The corresponding renormalized correlation functions involving N bulk, M surface fields and L\ surface operators 
\4> 2 {R,i) can be written as 

Gg ' M ' Ll) (p; m, u, v, c) = Z -("+m)/2 z -m/2 z ^ g{n , m , Li) (p . m ^ ^ VQ ^ y (31) 

In the present paper we concentrate our attention on correlation function G^°' 2,1 ^(p; m, u, v, c) involving two surface 
fields and a single surface operator insertion 2 (R;). 

It is well known [11] that the uv-singularitics of the correlation function q( n - m ,Li) can ^ e ac j sor bed through a mass 
shift 77Jq = m 2 +Sm 2 and surface-enhancement shift Co = c+5c. The renormalizations of the mass m, coupling constant 
u, v and the renormalization factor Z^ are defined by standart normalization conditions of the infinite- volume theory 
[27,26,28,29,2]. In order to adsorb uv singularities located in the vicinity of the surface, a surface-enhancement shift 5c 
is required. In this connection the new normalization condition should be introduced (see Appendix 1). Taking into 
account the normalization condition (A1.3) and expression for renormalized correlation function (3.1) it is possible to 
define the renormalization factor Za.i in the form 



, 7 aj^HOjmo.uo.vo.co)]- 1 
l z <P 2 J = z \\ 



dc 



(3.2) 

Co— Co (c,m,ii,v) 



It should be mentioned that renormalization factor = Z\Z§ is defined via the standard normalization condition 
(A1.2) (sec [11], [17]) 



= lim™£[G<°' 2 >( P )]- 1 . (3.3) 



Eq. (3.2) enables us considerably simplify the calculation of the correlation function G^ ' 2 ' 1 ) with surface operator 
</> 2 (R;) insertion. 

It should be noted that all Z factors in the d < 4 case have finite limits at A — > 00 (where A is a large-momentum 
cutoff). All factors mentioned above depend on the dimensionless variables u and v. Besides, the surface renormal- 
ization factors Z\ and Z$i depend on both u, v and the dimensionless ratio c/m. The last dependence on the ratio 
c/m plays the crucial role in the investigation of the crossover behavior from the special transition (c/m — > 0) to the 
ordinary transition (c/m — > 00). 

IV. EXPANSION OF THE CORRELATION FUNCTION NEAR THE MULTICRITICAL POINT 



As was indicated before, the main goal of the present work is to investigate the scaling critical behavior between 
special and ordinary transition and to calculate the crossover exponent <&. In this connection let consider the small 



deviations Ac — c — c* sp from the multicritical point. The power expansion of the bare correlation functions 
q(n,m)(j ) - mQj UQj Vq ^ Cq ) m terms of this small deviations Acq has a form 



G (Ar ' M) (p;r7i ,u , v ,c ) = 



Li=0 



(Ac ) j 



■GW M ' Ll )(p;m ,Mo^o,cl). 



(4.1) 



Based on Eq.(3.1), we rewrite the right-hand part of Eq.(4.1) in terms of the rcnormalized correlation functions 
and renormalized variable Ac = [Z^u, u)] _1 Aco and obtain 



Z -(N + M)/2 {ZirM/2G{N , M){p . ^ Cq) 

- £ ^4 Wl) (P; m,u,v). 



Li=0 



(4.2) 



The last equation in straightforward fashion define the correspondent renormalized correlation functions defined in 
the vicinity of the multicritical point 



G^ M \p;m,u,v,Ac) = ZT^^Z^G^Xvm^v^). 



(4.3) 



It is easy to see that these correlation functions depend on the dimensionless variable c = Ac/to. Thus, the correlation 
functions G^ ,M \p;m, u, v, Ac) satisfy correspondent Callan-Symanzik equations [31,11] 



m dm~ + v) lkL + ^ v ^Ih + — 2 — ^ v) 



M _ d 

+ Y??i P (w, u) - [1 + ??g(u, w)]c^r 



G^f' M) (p; m, u, u, Ac) = AGr, 



(4.4) 



where the inhomogcncous part AG^ should be negligible in the critical region similarly as that takes place in the 
case of infinite field theory. The resulting homogeneous equation differs from the standart bulk Calan-Symanzik (CS) 
equation [32-34] in that fashion it involves the additional surface related function rf^ and term — [1 + rj B (u,v)]c-^, 
where 



, , dlnZAu,v) n . .dlnZi(u,v) 
lnZ x (u, v) = [3 U (it, v) -^-^ + ft, (u, «) - 



FP 



i9u 



9u 



and 



rjc(u,v) = m 



dm 



dlnZj,2(u,v) dlnZj,2(u,v) 
lnZj,2(u,v) = f3 u (u,v) -2 \-/3 v (u,v) 



FP 



du 



dv 



FP 



FP 



(4.5) 



(4.6) 



It should be mentioned that functions /3 u (u,v), /3 v (u,v) and rj^u, v) appearing in (4.4) are the usual bulk RG 
functions. The symbol 'FP' indicates that the above value should be calculated at the infrared-stable random fixed 
point (FP) of the underlying bulk theory. 



V. SCALING CRITICAL BEHAVIOR AT THE MULTICRITICAL POINT 



The asymptotic scaling critical behavior of the correlation functions can be obtained through detailed analysis of 
the CS equations of Eq. (4.4), as was proposed in [32,35] and employed in the case of the semi-infinite systems 
in [20,11,21]. Our present investigations of the scaling critical behavior are in complete analogy with the scheme 
mentioned above [35,11] (see Appendix 2). Taking into account the scaling form of the renormalization factor Z^i of 
Eq. (A2.1) and the relation m ~ t" , we obtain for Ac and for the scaling variable c the next asymptotic dependences 

Ac-m-^'^Aco, Ac-t-^^'^Aco (5.1) 

and 

c~m- (1+ " E(u * x)) Ac , c~r"*Ac , (5.2) 



where 



$ = i/(l + »75(u*,t;*)) (5.3) 

is the surface crossover critical exponent. Eq. (5.2) explains the physical meaning of the surface crossover exponent as 
a value which characterises the measure of deviation from the multicritical point. The second equations in Eqs.(5.1) 
and (5.2) indicate about non-analitic temperature dependence of the renormalized surface-enhancement deviation Ac. 
Taking into account the above mentioned results from the CS equation we obtain the next asymptotic scaling form 
of the surface correlation function G^ ' 2 - 1 

sp 

G<°' 2 >(p; m , uo, v , c ) ~ m~^G^ 2) (^; 1, u\ v* , rrT^Aco) 

m 

~ T -^G(pt- v ;l,T-*Ac ), (5.4) 

where 7^ = v(l — rj\\), is the local surface susceptebility exponent and rjf, = rjl P + r/ is the surface correlation exponent 
[36] . It is easy to see, that the asymptotic scaling critical behavior of the surface correlation function for the systems 
with random quenched bulk disorder is characterized by the new crossover exponent $(it*, v*), which belongs to the 
universality class of the random model. In the next section, we will calculate the surface crossover exponent $ of the 
semi-infinite systems with random quenched bulk disorder. 



VI. THE PERTURBATION SERIES UP TO TWO-LOOPS. 



According to the Eqs.(5.3) and (4.6) the calculation of the crossover critical exponent <E> is connected with calculation 
of the renormalization factor via Eq.(3.2). The usual bulk ^singularities which are present in correlation function 
[G( 0,2 )(0)] _1 can be removed by the method similar to those repurted in Refs. [28,30,11,16] with help of standard 
mass renormalization procedure. 

The second step of our calculation is to remove the uv divergences which are connected with the presence of the 
surface in the system. The surface uv-singularities of the inverse surface correlation function [G^°' 2 ^(0)] _1 can be 
removed by performing the surface enhancement renormalization which is defined by Eq.(Al.l). For convenience we 
can rewrite the normalization condition of Eq. (Al.l) in the form 

Z ] \[G^\fr,m Q ,u ,v Q ,c Q ))- 1 =m + c. (6.1) 

for inverse unrenormalized surface correlation function [G(°' 2 )(0)] _1 . Performing the differentiation of the above men- 
tioned normalization condition with respect to and taking into account Eq.(3.2) we obtain for the renormalization 
factor Zm2 the next equation 

= % m 

where cq = c + 5c and 

5c = (Z' 1 - l)(m + c) + cr o (0; m, c = c + 5c). (6.3) 

Here oo(0;m, Co) denotes the sum of loop diagrams of all orders in [G^°' 2 ^(0; m, Uo, vq, Co)] -1 (see [11,17]). Among 
them o 1 corresponds to the one-loop graph, 02 denotes the melon-like two-loop diagrams 



G 
G 

°2 = V G 



(6.4) 



k 2 =0 



03 and 04 represent the reducible and irreducible two-loop diagrams in [G^ 0,2 ^ (0; m, Uo, Vq, Co)] -1 , respectively. Here 
the full lines with labels " G" denote the full free propagator of Eq. (2.4) . Eq. (6.2) can be resolved by using the method 
of sequential iteration. As a result of the first order of the perturbation theory (one-loop approximation) at general 
spatial dimensions d < 4, we obtain 

7 (l)_ 1 r X 7rj[f r(e) i+£ 3-e_l + e 3 + e 1 

Z ^ - 1 " Y 16 r(4i) [1 " 2 2Fl( ^~' ~2~> —> 2 )] ' (6 - 5) 



where 2-F\(...) is the hypergeometric function and coefficient T\ = !1 ^-Uo+v includes dimensionless effective expansion 
parameters u = u m~ e and v = v m~ e which are identified through the standart vertex renormalization. In the 
case of massive theory the value of e is not necessary restricted to be a small and the above expressions hold for any 
relevant dimensions d < 4. The presence of the Euler gamma function T(e) indicates the existence of dimensional 
poles 1/e when e — > 0. 

According to Eqs.(3.2) and (6.5) for r/ s at a one-loop order in the case of general dimensions up to d — 4 we obtain 

* " T i^"^im - 2*^(^1^; (6.6) 

At the random fixed point of order 0(y/e) [37](iCtu* = — 3^/||, K 4 v* — 4^/||, where the geometric normalization 
factor K 4 = 1/(8tt 2 ) ), Eq.(6.6) in the limit e -> leads to 

lim^ = lim„ c - = -y|. (6.7) 

This result coincides with that obtained by Ohno and Okabe [37]. At e — > for the surface crossover exponent $ in 
one-loop calculations we obtain 

$ = ---4/—. (6.6 
2 4 V 53 v 

In the case of three spatial dimensions (d = 3) the renormalization factors Z\ and Z$i are finite and their one-loop 
expressions do not coincide. At one-loop order, we obtain 

?7 e w -0.596 and $ = 0.286. (6.9) 

In the next order of the perturbation theory we restrict our attention only to the case of d = 3 dimensions. Thus after 
the surface-enhancement renormalization and performing the Feynman integrals in analogy with [11,17] and carrying 
out the vertex renormalizations of bare dimensionless parameters uo = uq/ 87tto and vo — v o /87rm 

W1 n + 8_ _, 

u = u(l H — u + v), 

6 

vo =v(l + -v + 2u), (6.10) 

we obtain a second-order series expansion for the renormalization factor Z^i in terms of new renormalized coupling 
constants u and v, 

Ztf(ii,v) = l + ^|-^(/n2- ^)u + (ln2- + ^^C{n)u 2 + 2C{n)uv + C{l)v 2 , (6.11) 

where C{n) is a function of the replica number n, defined by 

C(n) =A-B- \ln2 + ^ln 2 2 + ^±1, (6.12) 

and A = 0.202428, B = 0.678061 are integrals originating from the two-loop melon-like diagrams. Combining 
the renormalization factor Z^i with the one-loop pieces of the (3 functions (3 a (u 7 v) — —u(l — [(n + 8)/6]u — v) and 
Pv(u, v) = —v(l — § v — 2u) according to Eg. (4. 6), we obtain the desired series expansion for rj s , 

v) = -2^±|(M - \)u \{ln2 - \)v - 8[3^±^D(n)u 2 + 2 -^uv + (6.13) 

D{n) = A-B + 2±l ln * 2 _ H±l ln2 + 1!^, (6 . 14 ) 

and renormalized coupling constants u and v, normalized in a standard fashion u = [(n + 8)/6]u and v = |u. In 
common, Eq.(6.13) gives a result for the model with the effective Hamiltonian of the Landau-Ginzburg- Wilson type 



where 



with cubic anisotropy in the semi-infinite space (2.2) with general number n of order parameter components. Our 
calculations are connected with the investigation of the critical behavior of semi-infinite random Ising-like systems by 
taking the replica limit n — > 0. Hence, we obtain 

1, Is 2, L , 3 , N , D(0) D(l) 9l 

r,- c = --{ln2 - -)u - -{ln2 - -)v - S[^D{0)u 2 + -U-uv + -^« 2 ]. (6.15) 

The knowledge of r\c gives access to the calculation of the crossover critical exponents <& via the scaling relation of 
Eq.(5.3). Besides, we can calculate the critical exponents at\ and an of the layer and specific heats via the usual 
scaling relations [4] 

ai = a + u- 1 + $ = 1 - v{d- 2 -r} 5 ), an = a + v - 2 + 2$ = -v[d - 3 - 2jy s ]. (6.16) 

Above critical exponents should be calculated at the standard infrared-stable random fixed (FP) point of the 
underlying bulk theory [38] u* = —0.60509 and v* = 2.39631, as it is usually accepted in the massive field theory. 

VII. NUMERICAL RESULTS 

For each of the surface critical exponents mentioned above and crossover exponent $ we obtain from Eq.(6.15) at 
d = 3 a double series expansion in powers of u and v truncated at the second order [39]. In order to perform the 
analysis of these perturbative series expansions and to obtain accurate estimates of the surface critical exponents a 
powerful resummation procedure must be used. One of the simplest ways is to perform the double Pade-analysis [23] . 
This should work well when the series behave in lowest orders " in a convergent fashion" . Another way is to perform 
the double Pade-Borel analysis [24] for these series. The usage of the Pade-Borel resummation procedures are posiblc 
in the case when the terms in the series are alternating in sign [40]. The results of our calculations at Jug random 
fixed point [38] are represented in Table 1. The quantities O1/O2 and Ou/O^i represent the ratios of magnitudes 
of first-order and second-order perturbative corrections appearing in direct and inverse series expansions. The larger 
(absolute) value of these ratios indicate about the better apparent convergence of truncated series. 

The values [p/q] (where p, q = 0, 1) in Table 1 are simply Pade approximants which represent the partial sums of 
the direct and inverse series expansions up to the first and second order. The nearly diagonal two- variable rational 
approximants of the types [11/1] and [1/11] give at u = or v = the usual [1/1] Pade approximant [23]. The results 
of Pade-Borel analysis of the direct R and the inverse R~ x series expansions give numerical estimates of the surface 
critical exponents with a high degree of reliability. As it is easy to see, the most reliable estimate is obtained from 
the inverse series expansion for the surface critical exponent ct\, which represent the best convergence properties. 
Substituting this value a\ — 0.211 together with the standard bulk value v = 0.678 into the scaling relations (5.3) 
and (6.16), we have obtained the remaining critical exponents that are present in the last column of Table 1. The 
deviations of these estimates from the other estimates of the table give a rough measure of the achieved numerical 
accuracy. 

The results of the similar analysis of the perturbative series expansions of the surface critical exponent at random 
fixed point u* = —0.6524, v* = 2.4203 [41] are presented in Table 2 for comparison. In a similar way the most reliable 
estimate is obtained for inverse series expansion of a\. The results of substituting of a\ — 0.208 and v — 0.679 into 
scaling laws (5.3) and (6.16) are presented in the last column of Table 2. As it easy to see from comparison of the 
results Table 1 and Table 2, the difference in the ways of the (3 functions resummation have not essential influence on 
the values of the surface critical exponents. The difference between final results of Table 1 and Table 2 are 1.2% for 
n s , 1.4% for cm, 1.8% for a n and 0.2% for $. 

For evaluation of the reliability of the results obtained in the two-loop approximation, we have performed additional 
calculation of the surface critical exponents from ot\ = 0.211 and six-loop perturbation theory results [2] for bulk 
critical exponent of the correlation length v = 0.678(10). We have obtained rj e = —0.164, an = —0.222 and for 
surface crossover critical exponent $ = 0.567. This indicates about good stability of our results obtained in the 
frames of two-loop approximation scheme. 

VIII. SUMMARY 

We have studied the crosover critical behavior between special and ordinary surface transitions of three-dimensional 
quenched random semi-infinite Ising-like systems. We find that the asymptotic scaling critical behavior of the surface 
correlation function for the systems with random quenched bulk disorder is characterized by the new crossover 



exponent &(u*,v*), which belong to the universality class of the random model. We have calculated the crossover 
critical exponent $ and critical exponents of the layer, ai, and local specific heats exponent, an, by applying the field 
theoretic approach directly in three dimentions up to the two-loop approximation. We performed a rational double 
Pade and double Pade-Borel analysis of the resulting perturbation series expansions for the surface critical exponents 
in order to find the best numerical values. The final numerical values of the surface critical exponents a\, an and 
crossover exponent $ for the systems with quenched random bulk disorder are 

on =0.211 ±0.003 an = -0.222 ±0.004 $ = 0.567 ± 0.001. (8.1) 

These values evidently different from their counterparts of pure Ising systems [30,11] 

ai = 0.279 an = -0.182 $ = 0,539. (8.2) 

We suggest that the obtained results could stimulate further experimental and numerical investigations of the 
surface critical behavior of random systems. 
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APPENDIX 1 

In order to specify 8c, Z\ and Z^, we require that [30,11] 

G^' 2 \p;m,u,v,c) = 1 , (Al.l) 
P =o m + c 



dGft ' 2 ' (p; m, u, v, c) 



dp 2 



p=0 



2ra(m + c) s 



(A1.2) 



and correspondent normalization condition for the correlation function G^ 0,2,1 ' with the insertion of the surface 



operator ^<f) 2 



G , l ' 2 ' 1 '(p;m,M,«,c) 



1 



p =o (m + c) 2 



(A1.3) 



Eq.(A1.3) is motivated by the fact that the bare correlation function G^ - 2 ' 1 ) (0; m , u , v , c ) may be written as 
derivative — ^-G^ 0,2 ^(0; m , uq, vq, cq). This equation simplifies considerably the calculation of the correlation function 



1 A2 



with insertions of surface operator 2 

From Eq.(Al.l), it is easy to see that the special point is located at m = c = 0, because at this point the divergence of 
the bulk and the surface correlation length and susceptibility is observed. At c = the surface normalization conditions 
are simplified and yield Co = c* p . This point corresponds to the multicritical point (wg c , c* p ) at which special transition 
takes place. On the other hand, the above mentioned equation implies also that the surface correlation length and 
the susceptibility are finite at the ordinary transition, because in this case we have c > when m — > 0. This latter 
case corresponds to the situation when the surface remains " noncritical" at the bulk transition temperature. 



APPENDIX 2 



As it is usually accepted in the massive field theory, the variable m is identified with the inverse bulk correlation 
length and is proportional to r v , where r = (T — T c )/T c . Following the scheme proposed in [35], we can perform 
the integration of Eq.(4.5), Eq.(4.6) and expressions for the RG functions (3 u (u,v), (3 v (u,v) and f]<j,(u, v). This gives 
the following asymptotic dependencies at m — > 



= P U {U*,V*), 

u v = p' v (u*,v*), 

(A2.1) 

As follows from these expressions, the variables u and v deviate from their fixed values u* and v* by different 
scaling laws with various values of lo u and u) v . The scaling laws (see Eq.(A2.1)) have the similar form as in the case of 
the pure systems [11], but renormalization factors are characterized by another values of the critical exponents which 
belong to the universality class of random model. 

TABLE I. Surface critical exponents involving the RG function rjc at the Jug fixed point u* = —0.60509, v* = 2.39631 
(two-loop order) 



\u — it*|^m"", where 
\ v — v *\ <~ m w ", where 
Z<j, ~ m n , where 



exp 


o, 


(hi 

o„ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


[R] 


K 1 


/(ai,">»7) 


Ve 


-0.8 


-1.5 


0.00 


-0.574 


-0.365 


0.150 


-0.152 


-0.281 


-0.268 


-0.313 


-0.280 


-0.164 


Oil 


-1.7 


-7.6 


0.50 


0.051 


0.190 


0.312 


0.220 


0.201 


0.213 


0.185 


0.211 


0.211 


an 


-1.1 


-2.8 


0.00 


-0.574 


-0.365 


-0.036 


-0.268 


-0.324 


-0.306 


-0.351 


-0.313 


-0.222 


$ 


-0.5 


-0.5 


0.5 


0.375 


0.389 


0.652 


0.658 


0.451 


0.452 


0.444 


0.445 


0.567 



TABLE II. Surface critical exponents involving the RG function rjc at the fixed point u* = —0.6524, v* = 2.4203 (two-loop 
order) 



exp 


Oi 

o, 


o«. 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


[R] 








-0.82 


-1.55 


0.0 


-0.570 


-0.363 


0.124 


-0.168 


-0.287 


-0.272 


-0.318 


-0.283 


-0.166 


Cil 


-1.82 


-9.56 


0.5 


0.054 


0.191 


0.300 


0.215 


0.197 


0.209 


0.182 


0.208 


0.208 


an 


-1.12 


-3.08 


0.0 


-0.570 


-0.363 


-0.060 


-0.278 


-0.331 


-0.311 


-0.357 


-0.317 


-0.226 


$ 


-0.47 


-0.5 


0.5 


0.376 


0.389 


0.641 


0.643 


0.449 


0.450 


0.442 


0.444 


0.566 
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